Charge order in LuFe204: antiferroelectric ground state and coupling to magnetism 
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X-ray scattering by multiferroic LuFe204 is reported. Below 320 K, superstructure reflections in- 
dicate an incommensurate charge order with propagation close to (|§|)- The corresponding charge 
configuration, also found by electronic structure calculations as most stable, contains polar Fe/O 
double-layers with aniiferroelectric stacking. Diffuse scattering at 360 K, with (|-|0) propagation, 
indicates ferroelectric short-range correlations between neighboring double-layers. The temperature 
dependence of the incommensuration indicates that charge order and magnetism are coupled. 

PACS numbers: 71.30.-|-h, 77.80.-e, 75.80.+q, 61.05. C- 
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Materials where ferroelectricity or dielectric behavior 
is coupled to magnetism have the potential for novel ap- 
plications and presently receive a lot of attention [l| . A 
new type of ferroelectricity, originating from charge or- 
der (CO) and seemingly coupled to magnetism, has been 
proposed to occur in LuFe2 04 containing triangular Fe/O 
enerating a lot of interest in this ma- 
Ferroelectricity is thought to 
arise from a particular arrangement of Fe^"*" and Fe'^+ 
within the Fe/O double-layers, making these intrinsically 
polar. However, different reported [8|, |9[ superstructure 
reflections indicate an incomplete understanding of the 
CO, and the full three-dimensional (3D) charge configu- 
ration has yet to be established. The latter determines 
the stacking of the polarizations of the individual double- 
layers, and thus the net polarization of the material. 

Further, although magnetism in LuFe204 occurs in the 
CO state to date there is no direct observation of coupling 
between CO and magnetism to date. Such an observa- 
tion would be important in establishing the mechanism 
by which ferroelectricity, dielectric behavior, and mag- 
netism couple via the underlying CO. We have recently 
grown LuFe204 crystals with magnetic transitions of un- 
precedented sharpness, on which neutron scatte ring al- 
lowed the first refinement of a 3D spin structure [lO| . 

Here, we present a study of the CO superstructure 
by synchrotron x-ray scattering. We propose a com- 
mensurate approximation for three domains of the in- 
commensurate CO configuration, with propagations (in 
hexagonal notation) close to the symmetry-equivalent 
directions (Hf), (iM)i and (||§), corresponding to 



an anijferroelectric [11[ ground state that we also iden- 
tify by first-principles density-functional theory (DFT) 
as having the lowest energy. In contrast, short-range 
charge correlations above the CO temperature were 
found to correspond to a ferroelectric CO configuration 



with (-I^O) propagation. Further, we provide evidence 
for a coupling between CO and magnetism involving pri- 
marily the incommensuration of the CO. Our results un- 
derline the importance of near degeneracy of the CO in 
LuFe204 and provide an essential microscopic basis for 
magnetoelectric coupling previously proposed [3,l3|, 4, 5]. 

We studied single crystals from the same batch as those 
in [10|. X-ray scattering was performed on the (001) sur- 
face of a crystal with a mosaic of 0.02(1)° at beam line 
6IDi? of the Advanced Photon Source, using 16.2 keV 
photons. All scattered intensities are normalized to an 
ion chamber monitor. Structural Bragg reflections, de- 
noted s, are consistent with the reported [I4I rhombo- 
hedral _R3?Ti structure. Specific heat was measured with 
commercial equipment. 

Superstructure reflections (Fig. [T]) appear below a 
sharp feature in the specific heat (inset) at Tco ^ 320 K 
indicating the CO transition [13]]. We found two sets of 
strong superstructure reflections at the same hki posi- 
tions as reported in [9|, one set near (ggf) with o an 
odd integer, and satellites to (00^) at ±{ttO) and (not 
previously reported) two symmetry-equivalent directions 
[r'^0.03, see Fig.[2ja)]. As proposed in [o!], the positions 
of these refiections are consistent with three CO domains 
(labelled A, B, and C hereafter) with the symmetry- 
equivalent propagation vectors p^ = (i + (5,i + 5,|). 
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5 ^ 0.003 is very small, but 5 7^ implies incommensu- 
rate CO. The reflections near (^^f ) correspond to s±p., 
and several additional, much weaker, types of reflections 
visible in Fig.[l]have positions consistent with s±np , to 
be discussed in detail elsewhere. The reflections shown 
in Fig. [2] may be considered as higher harmonics as well, 
likely s±9p. (t = 9(5), and are in any case a consequence 
of a bimodal charge distribution of the CO configuration 
obtained below. 
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FIG. 1: (Color online) Scattered intensity at (hhl) at 130, 
200, 295, and 360 K. Lowest panel: {^^£) cut at 360 K with 
fit by Gaussians. Strong structural reflections reach inten- 
sities > 10 counts/monitor, strong superstructure reflections 
~40. The higher overall intensity at 295 K is due to wider de- 
tector slits being used. Horizontal lines mark hh = ^. Insets: 
Specific heat near 320 K; sketch of hki positions of maxima in 
diffuse scattering at 360 K, with symbols indicating domains 
as in Fig. EJa) . 



To estimate domain populations we collected x-ray 
diffraction data on a second crystal 14], using a CuKa 
diffractometer [Fig. [3fa,b)]. Intensities of superstructure 
reflections close to (Hf) and equivalent directions are 
shown in Fig.[2Ia). Reflections at positions s±p^, sips, 
and s ± pc are indicated by ■, ▲, and T, respectively. 
Taking into account that R3m structural reflections oc- 
cur only for —h + k + l = 3n with n integer, each observed 
superstructure reflection is associated with one particular 
domain. All measurements are consistent with domain A 
being almost unpopulated, and domain B roughly twice 
more populated than domain C. The intensity ratios of 
satellite pairs around (0,0,19.5) [Fig. [3I^b)], each pair 
arising from one domain, indicate the same domain popu- 
lations, supporting the hypothesis that they are resulting 
from the same CO configuration. The consistent domain 
populations obtained from the intensities confirms the 
PA, Pb, PC, propagation type proposed above. 

To determine the CO configurations corresponding to 
propagations P a, Pb, and pc, we performed represen- 
tation analysis |l5l. Ilq . Similar to the spin order deter- 
mined earlier ilOl, there are two allowed irreducible CO 




FIG. 2: (Color online) Incommensurate satellites, (a) In- 
tensity at (/ifcl3.5) at 130 K (the two strongest peaks reach 
~19counts/mon). (b) Detail around (0,0,13.5). (c) Intensity 
around {^114.5) (marked by o). 



representations i) and ii), corresponding to either same 
or different valence for the 1 and 2 Fe sites [Fig. [2Ic,d)] 
of the primitive cell. For simplicity, our following discus- 
sion assumes a commensurate approximation (5, r — > 0) 
[18| in which the propagation vectors become ( H | ) and 
symmetry-equivalent wavevectors. Both representations 
contain sites with two different magnitudes of the Fe va- 
lence difference from average, in clear contrast to the re- 
sults of Mossbauer spectroscopy, which imply a bimodal 
valence distribution |17j. A bimodal charge distribution 
is only obtained by adding a uniquely defined (00 1) rep- 
resentation. Case i) leads to a CO configuration in which 
the Fe/0 double-layers are not charge neutral. This is 
physically very unlikely, given the separation of neigh- 
boring double-layers by ~ 6 A. Case ii) leads to a single 
feasible configuration with overall neutral double-layers, 
see Fig. ^c). Shown is the configuration for domain A, 
accommodated in a V^ x V3 x 2 supercell. Note that 
the CO lowers [19| the crystallographic space group sym- 
metry to monoclinic C2/m. For this CO configuration, 
structure factor calculations *20| confirm the observed re- 
flection pattern, including both sip. and s±(00|). The 
latter, in the commensurate approximation, correspond 
to the reflections shown in Fig. [UJa), again conflrming 
the association of both these types of reflections with the 
same CO. As a further check of the proposed CO [Fig. 
[He)], we used this solution to re-analyze the 220 K neu- 
tron data presented in [10|. Combining the CO model 
with the previous model of spin order, magnetic scatter- 
ing at both i^^n) and (|-|§) could be refined successfully 

For the above CO model, the configuration in each 
individual double-layer is polar as proposed in [3]- How- 
ever, the stacking of the polarization of the six double 
layers of the supercell [Fig. [3{c)] is anizferroelectric with 
no net polarization. To confirm this antiferroelectric con- 
figuration, we carried out DFT calculations for the ferri- 
electric and antiferroelectric CO within a \/3x-\/3x2 cell 
(in addition to the y/^x^/Sxl calculations of [(|), that in- 
deed show that the antiferroelectric CO configuration is 
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FIG. 3: (Color online) Commensurate approximation of the 
charge configuration, (a) Integrated intensities of (ii|) type 
reflections. Reflections associated with different domains are 
indicated by different symbols (see text), (b) scattered inten- 
sity at (/ifef9.5). (c,d) Charge configuration for domain A for 
T < Too (c) and T > Tco (d, short-range). Shown are the 
Fe ions (see legend) of the 6 (3) Fe/0 double-layers of the 
\/3x\/3x2 (xl) supercell. The direction of the local polariza- 
tion for each double layer is indicated by an arrow. The two 
sites of the primitive cell are numbered (see text). 



more stable than the ferrielectric one [by 3.2 meV/fu (for- 
mula unit)], and hence corresponds to the ground state. 

Our result raises the question how the remanent polar- 
ization indicated by pyroelectric current measurements 
[2| could be explained. A simple explanation would be 
that the samples are different, since the strong oxygen 
stoichiometry-dependence of physical properties is well- 
known. However, this is unlikely to apply, for two rea- 
sons. First, superstructure reflections published in the 
same paper as the polarization results [2] and other pa- 
pers by the same group (e.g., |3]), are consistent with 
(ii|) propagation. Our structure factor calculations in- 
dicate that all CO configurations with non-zero net polar- 
ization have {■^hn) reflections of the same order of mag- 
nitude as (^jf ) reflections, whereas no (j^n) reflections 
are reported in [2, ISj. Second, preliminary pyroelectric 
current measurements with a similar protocol as in fZ| 
suggest a similar remanent polarization as in 2] . 

Before presenting an alternative explanation, we turn 
to the diffuse scattering observed above Tco- Heat- 
ing through Tco, the (Hf) and all satellite reflections 
weaken [Fig. [4l[c,d)], as expected, but {■^^n) reflections 
gain in intensity [Fig. Hljb)]. At 360 K there is still con- 
siderable diffuse scattering present around the (^|^) rod, 
see Fig. [TJ A cut in hM along {■hh^^) reveals broad max- 



ima now at integer £ positions, well described by strongly 
overlapping Gaussians (lowest panel). The positions of 
these diffuse peaks are sketched in an inset to Fig. [T] 
They are consistent with (^ — 5',^ — i5',0) propagation 
and equivalent. The width of the Gaussians in Fig. [T] in- 
dicates that CO correlations in c direction (^c) extend to 
2—3 double-layers [23|, i.e. they are 3D rather than 2D, 



in contrast to [Qj]. Representation analysis for the above 
propagation again yielded two irreducible representations 
[23(, one rejected as unphysical due to non-charge-neutral 
double-layers. The remaining representation in commen- 
surate approximation leads to a CO conflguration [Fig. 
|31[d)] with v^x-s/Sxl cell and ferroelectric stacking of the 
polarization of the double-layers. Our DFT calculations 
using the \/3x y/Sx 2 cell show that this ferroelectric con- 
figuration is less stable than the antiferroelectric config- 
uration (i.e., the ground state), but only by 13.4meV/fu 
or 3% of the overall CO gain. 

Although with £,c S ^ the description in terms of CO 
configurations may seem somewhat questionable, the re- 
sult suggests that the high T correlations favor a ferro- 
electric arrangement between neighboring double-layers. 
It is thus surprising that, upon long-range charge or- 
dering, the configuration established is not ferroelectric. 
Given the small energy-differences between antiferro— 
and ferroelectric configurations it seems likely that cool- 
ing with an electric field applied may stabilize long-range 
order of the ferroelectric CO configuration, which would 
explain the remanent polarization observed (only) after 
cooling in an electric field. This idea should be tested by 
scattering experiments with electric fields applied in-situ. 

In zero electric field, the ground-state has the same 
basic CO configuration with no net polarization at all T 
below TcOy but specific features vary with T. Here, we 
focus on the temperature dependence of the incommen- 
suration of the CO. As in Fe20B03 [2J], the incommen- 
suration r [25| changes with T [Fig.HJe)], but in LuFe204 
it is much smaller and present at all T. In the fluctuation 
regime, r decreases rapidly upon cooling below Tco, sug- 
gesting that the incommensuration may be an important 
factor in stabilizing either ferro- or antiferroelectric CO. 
This trend is reversed at Tjv ~ 240 K, where t ^ 0.025 has 
a minimum: The appearance of the seemingly commen- 
surate [10| magnetic order results in forcing the CO to 
become more incommensurate. In the center of the satel- 
lite rings [Fig.[2Ua,b)], weak additional reflections develop 
upon cooling through Tn [Fig. Hlja)]. The emergence of 
these additional reflections is not completely understood 
at present, but provides an additional indication of cou- 
pling between the CO and the magnetism. 

At the low T transition Tj^t ~ 175K reported in [10| 
the width and intensity of these additional reflections 
change, and a moderate impact on various other reflec- 
tions and t{T) is readily visible in Fig. H) While a de- 
tailed discussion of the impact of the low T transition on 
CO and other degrees of freedom is beyond the scope of 
this work, we like to point out that i) the basic CO con- 
figuration remains the same upon cooling through Tlt, 
and ii) the CO correlation along c is slightly (~ 7%) 
improved compared to 200 K, suggesting a better estab- 
lished CO. In contrast, the magnetic order is less estab- 
lished below Tlt as indicated by the broadness of various 
magnetic refiections and diffuse magnetic component re- 





FIG. 4: (Color online) T dependence of CO on warming, (a) (/i/il3.5) scans, (b) Scans ||(120) through (|il4). (c) (||13.5) 
rocking curves, (d) Peak intensities of three reflections vs T. (e) Incommensuration t{T). 



ported in [lO]. Thus the "final compromise" between CO 
and magnetism, established below Tlt on approaching 
the ground state, seems more favorable to the CO. 

In conclusion, we show by scattering experiments and 
DFT calculations that the charge-ordered polar Fe/0 
double-layers of LuFe2 04 have antiferroelectric stack- 
ing in the ground state. Upon long-range ordering, the 
high-T ferroelectric short-range correlations revert to this 
stacking. The incommensurate nature of the resulting 
CO is likely relevant for stabilizing either ferro- or antifer- 
roelectric charge configurations, and reveals the coupling 
of the CO to the magnetism, a coupling likely related to 
the observed [3| large magneto-dielectric effect. 
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